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Abstract 

We compute the boundary scattering amplitudes of the breathers of the super- 
symmetric sine-Gordon model using the fusion of the soliton-antisoliton pair scat- 
tering with the boundary with a known result of the soliton boundary scattering 
amplitudes. We also solve the boundary Yang-Baxter equation of the eight-vertex 
free fermion models to find the boundary reflection matrices. The former result is 
confirmed by the latter since the bulk S'-matrices of the breathers can be identified 
with the trigonometric limit of the Boltzmann weights of the free fermion models. 
Our dual approach can answer a few quesions on the relationships between the free 
parameters in the boundary potential and those in the scattering amplitudes. 
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1 Introduction and Motivation 



The study of the two-dimensional integrable models of quantum field theories and 
statistical models based on the Yang-Baxter equation (YBE) has provided important 
understandings of nonperturbative aspects of these models and technical tools for appli- 
cations to real physical problems. The YBE plays essential roles in establishing the inte- 
gr ability and solving the models. In the field theories, the YBE provides a consistency 
condition for the two-body scattering amplitudes (S-matrices) in the multi-particle 
scattering processes since the scattering is factorizable. With unitarity and crossing 
symmetry, the YBE can determine the S'-matrix completely. In addition, the correla- 
tion functions may be obtained by computing multi-particle form-factors. The lattice 
models which are defined by the Boltzmann weights can have well-defined transfer ma- 
trices if they satisfy the YBE and can be diagonalized by independent technologies, 
such as the algebraic Bethe ansatz to extract exact properties of the model. 

Recently there has been a lot of efforts in extending these approaches to models 
with boundaries. They are motivated by the fact that these models with the boundary 
have more applicability to real physical systems than those without one. For example, 
three-dimensional spherically symmetric physical systems can be effectively described 
on the half-line if s-wave element becomes dominant. One-channel Kondo problem, 
monopole-catalyzed proton decay are frequently cited examples. Also one can generalize 
the conventional periodic boundary condition of the statistical models to other types 
like the fixed and free conditions. 

The existence of the boundary adds new quantities like boundary scattering am- 
plitudes and Boltzmann weights, and one needs to extend the YBE to include these 
objects. The boundary Yang-Baxter equation (BYBE) (also known as the reflection 
equation) ]l|] plays the role of the YBE for the integrable statistical models [Q, ||] and 
quantum field theories Q| in the presence of a boundary; it is the necessary condition 
for the integr ability of these models. 

The super symmetric sine-Gordon model (SSG) preserves the integrability in the 
presence of the boundary || and the scattering amplitude of the SSG solitons with it 
has been computed ||. In this paper we compute the boundary scattering amplitudes 
of the SSG breathers, the bound states of the soliton and antisoliton pair, in two 
independent ways. The first one is to consider the the eight-vertex free fermion models 
with the boundary and to solve the BYBE. This is related to the SSG model since 
the breather S'-matrices are trigonometric limit of the free fermion models in certain 
regimes. The second is to compute the two-particle (a SSG soliton and an antisoliton) 
boundary scattering amplitudes and to take a limit of two rapidities so that they can 
form a bound state. This 'fusion procedure' can give independent check of our results. 
In addition, these two approaches can answer a few questions raised in the previous 
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study; how the super symmetry of the SSG lagrangian with the boundary can be realized 
in the scattering matrix context and how the parameters in the lagrangian and those 
in the boundary scattering amplitudes are related. 



2 SSG Breathers on a half line 

The action of the SSG model is given by 



S 



I 



dxdt 



-(<9 M 0) 2 -iip flip - -jp cos2 ^ ~~ 2m(cos 



(2.1) 



where is a real scalar field and ip is a Majorana fermion. [3 is a coupling constant and 
m is the mass parameter denoting the deviation from the massless theory. This theory 
has soliton spectrum \K^ b {Q)) where 'a&' and '±' are the RSOS spins (a, b = 0,^,1) 
and the topological charges ('+' for the soliton and '— ' for the antisoliton) , respectively 
and is the rapidity. The exact S-matrix of the SSG (anti-)solitons has the factorized 
form of ]7j 

Sssg(0) = s£§ g (0)®S S g(0). 

The first S'-matrix factor which acts on the supersymmetry (SUSY) charges is the 
RSOS S'-matrix for the tricritical Ising model perturbed by the $13 operator; 



Q 



(4) 
RSG 



dc 



(9) = U(9) (x$)~ L \y[x* 



' sinh [ - ) Sdb + sinh 



ITT 



, (2-2) 



for IK^O) + \K ab (9 2 )) -> \K dc (6 2 )) + |K c6 (^)) where X$ 



q- number [n] = (q n — q n )/(q — q x ) and q 



-i7r/4 



[2a+l][2c+l] ] 
[2d+l][2b+l] } WltJl 



(Fig.l). 




Fig. 1 Bulk S-matrix of the kinks 
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The second one is formally the sine-Gordon (SG) (anti)soliton S-matrix with the 
parameter given by 7 = 4/3 2 /(l — /3 2 /4-7r).[] The factorized form of the S'-matrix implies 
that the SSG soliton can be formally written as |K^(0)) = |K a6 (0)) (g> | ± (0)). The 
charge conjugation of the SSG solitons is defined by 

C\Kf b ) = |K£>. (2.3) 

For n < 87r/7 < n + 1, the second factor, the SG S'-matrix, has n poles in the physical 
strip corresponding to the SSG breathers. The threshold value of the SSG (3 to have 
any bound state is f3 2 = 4tt/3 compared with that of the SG, (3 2 = An. 

The bulk S-matrices of the breathers have been obtained by considering the residues 
of two solitons and two antisolitons scattering and taking appropriate limits on the 
rapidities. Due to the factorized form of the soliton S-matrix, the SSG breather S- 
matrices are also made up of two factors. The factors coming from the SG sector 
have been computed in || and they are completely diagonal since the masses of the 
SG breathers are non-degenerate. The second one comes from the four kinks scatter- 
ing fusion processes of the RSG(4) ||. It is obvious that this factor is non-diagonal 
since the SUSY makes the mass spectrum degenerate and the breathers form N = 1 
supermultiplets . 

Since two kinks can scatter only when they share a common RSOS spin, the two- 
kink states which form the SSG breathers can be written as 

\<m = ^(|K o i(0i)Ki 1 (0 2 ))-|K 1 i(0 1 )K |o (0 2 )) 

\<t>M) = -^(|K o i(0 1 )Ki o (0 2 )) + |K 1 i(0 1 )K |l (0 2 )} 

WniO)) = ^|(|K |o (0 1 )K o| (0 2 ))-|K |l (0 1 )K 1 x(0 2 )) 

l^nW) = 7 ^(|Ki o (0 1 )K o i(0 2 )) + |K |l (0 1 )K 1 i(0 2 )) 



(2.4) 



where the rapidities are related as 



= \ (0i + 2 ) , 0i - 2 = A6 



and 



a n = 4 / tan , A0 n = in — 



1 Notice that the S'-matrix of the SG model depends on 7 = f3 2 /(l — f3 2 /8iv). 
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with A9 n corresponding to the mass pole of the breathers. 

Notice that the bound states come in two sets distinguished by their superscripts 
1, 2 and only particles with the same superscripts can scatter. Each pair (ip^, <Pn) forms 
N = 1 supermultiplet (see Eq. flP]) ) and has the same bulk 5-matrices, |J 



/ sin + | sinh ( 



S(9) = p{6) 



\ 



-i sin ^ sinh | 



— | sinh 9 
sin cosh | 



sin cosh | 
— ~ sinh# 



i sin ^ sinh | \ 



sin yg- — § sinh ( 



/ 

(2.5) 

where the column and row are arranged in the order of ipffln, ^n^ni ^n^ni ^nfin f° r 
a = 1,2. The states in Eq.(2.4) are invariant under charge conjugation C, implying 
that they are real scalar particles and Majorana fermions. The least massive bound 
states (n = 1) are identified with the tp and tfi fields in the Lagrangian Eq.( |2.lD and, 
indeed, the above 5-matrix is identical to that obtained in [1C] if we identify sin ^ 
with /. Since there is only one fundamental field pair (ijj,(f)) in the lagrangian, the 
two sets of the bound states which have the same S'-matrices should be indentified, i.e. 
|V*(*)) = |V£(*)> and \<j>l{0)) |#(0)>. 

The function p{9) satisfies the unitarity and crossing relations 



p{9)p{-9) (sin 2 $ + S i n h 2 §) cosh 
p(9) = p(iir - 9) 



2 e 



(2.6) 



The minimum solution to these equations have been given as |1C] 

2i 



sinh 



-Z{Q)Z(m-9) 



(2.7) 



where 
Z{9) 



T(-i6/2ir) 
T(l/2 - «0/2tt) 



n 

i=i 



r(n7/16vr - i9/2ir + ^T^-nj/Wn - i9/2n + 1-1) 
r(n7/16vr - i9/2n + I + l/2)r(-n7/16vr - i9/2*rr + I - 1/2) 



T 2 (-i9/2Tr + l- 1/2) 

x r 2 (-^/27r + /- 1) 

Now we introduce a boundary potential which preserves the integrability. The SSG 
boundary potential that gives conserved charges at the first order has been derived as 



B(4>, ij), if}) = A cos 



(2i 
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With the assumption of the complete integrability, one can use the BYBE to solve this 
model. For the purpose, we use an important property of the S'-matrix in Eq.(|2.5|) that 
it satisfies the free fermion condition Jy], [|. Therefore it should be a special limit of 
the Boltzmann weights of the eight-vertex free fermion model given in the Appendix. 
Indeed, consider the regime \h\ > 1 and take the following trigonometric limit: 

1 

k — > with k cosh 5 = k sinh 5 = - 



sin^ 



We obtain the breather S-matrix if u = —id. An immediate consequence of this observa- 
tion is that the boundary S'-matrix of the SSG breathers are given by the trigonometric 
limit of the boundary Boltzmann weights derived in the Appendix. The trigonometric 
limit of Eq.( frT7 ) becomes 




cosh§G + (#) -zsinh§G_(0) -iesinhtf 

-isinh# coshfG+(0) + isinh §G_(#) 

ea± + a T 2 



(2.9) 



G±(6) =a±- 2 »7 — sinh ^ ( 2 - 10 ) 



and 



sin — e 2 



The overall factor p(6),1Z(0) are functions that ensure the unitarity and crossing sym- 
metry of the S- and R- matrices, respectively. 

The equations that determine 1Z{0) are given by 



n{6)n{-6) (cosh 2 \G\ (9) + sinh 2 §G 2 (6») + e sinh 2 0) = 1 
p(26)TZ(^ + 0) cosh (sin ^ - ie sinh 0) = ft(f - 0) . 
As usual we define 

= ft o (0)fti(0) 



(2.12) 



such that 



where 



fti(0) =Ki(iir-e) 

Tli(0)Ki(-9) (cq + ci sinh 2 f + c 2 sinh 4 f 



(2.13) 



2 e (a + + ea_) 2 i n (a + + ea_ 



cq = a, , ci = ^-^ h 2e , c 2 



sm ^ ( e - sm ^) 
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and 

p(29)K { l f + 6) cosh 9 (sin ^ - iesinh6>) = ft (f - 0) 

(2.14) 

TZ o (9)n o (-9)cosh0 = 1. 

The factor TZ\(9) carries information about the boundary conditions that are deter- 
mined by the free parameters a±, and its minimum solution is given by 

n 1 (9) = —a( X ,9)a(r ] ,e) (2.15) 

where the function a(x, 9) is an infinite product of T function defined as 

n( X , vr/2 + i9)U(- X , vr/2 + i9)U( X , -*/2 - i9)U( X , -*/2 - id) 



n(x,-*0) = n 



n 2 ( X ,7r/2)n 2 (- X ,vr/2) 
T{1 + x/vr + i9/2ir) 



l= \ F(l + 1 + x/vr + i0/2ir) ' 

with the parameters %, r\ defined by 

cos -2 x + cos -2 rj = ci/co , cos~ 2 %cos~ 2 r] = C2/C0 ■ 

The relations that determine TZq(9) can, similarly, be solved with minimum solution 
given by 

= gU (2 . 16) 

where 

r(l - l + en7/36vr + 1/4 + iO/2ir)T(l - en-f/36ir - 1/4 - ^/2tt) 



m = n 



L x r(-Z + en 7 /367r + 1/4 - i9/2ir)T(l + 1 - en 7 /367r - 1/4 + z(9/2tt) 
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3 Soliton Fusion 



In , we studied the scattering theory of the SSG model on a half line based on its 
soliton states. Essentially, in the presence of a boundary, integrability of the SSG model 
requires that the boundary S-matrix of the solitons satisfy the BYBE, Eq.( |5.2D . We 
can assume naturally that this boundary S'-matrix is also factorized into two parts: 

Rssg(0) = R { il G (0) ® Rsg(9) (3.1) 

where i?sc and Rrsq are the SG and RSOS(4) boundary scattering matrices, respec- 
tively. Writing it in this form, we can solve the BYBE separately. The SG part has 
been obtained in [Q, [l2j and the kinks part has been found in @ to be of the form 
(Fig.2) 

R a bc (9) = K{9) (x£)~& [5 b ^ c (e) + 5 bc (s b _ l/2>a U a (9) + S b+1/2ja D a (9))] (3.2) 
where o, b, c are the RSOS(4) spins. 




Fig.2 Boundary kink S'-matrix 



The explicit solutions are given by 

X& = S X{ , ^^ = -4r + ^cosh|, _^^ = _# -Ccoshf , 
Y2 sinh- vi smh| z 

A Q1 2 A Q1 2 



Di(fl) _ I- A sinh | 
W) ~ l + ^sinh| 



(3.3) 



with A,B,C being the free parameters of the boundary, and the off-diagonal terms 

i i 

X^jX^q are independent of the spectral parameter and differ from each other by a 
gauge factor s. The overall function 71(9) that guarantees boundary crossing and 
unitarity is given in || . 
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As we have shown that the bound states of the RSG(4) kinks give rise to the 
ip n {Q), fields, whose bulk scattering matrix is given by Eq. fl2.5Q , the scattering of 

these ip n , <j) n fields with the boundary is governed by the matrix given in Eq.( |2.9| ). Hence 
we expect that by fusing the boundary scattering matrices of the kinks given above, we 
should reproduce Eq.(2.9). Before proceeding with the computation, it is worth recall- 
ing that, the fermionic and bosonic bound states come in two types; (-0^, 4>n) , (V'ni 0n)i 
which are, nevertheless, identified in the bulk since the scattering matrices of these 
two sets of particles have exactly the same form and are hence indistinguishable from 
each other. It would be natural to wonder whether the same holds true in the pres- 
ence of a boundary. In fact, more intriguingly, notice that the boundary scattering 
matrices of the kinks carry more than one free parameters, while the scattering matrix 
given in Eq.( |2.9| ) has only one. From the fusion equation, it is clear that the .R-matrix 

of (?Pni $n) ■> which will be built out of Uq,D\, will contain one free parameter as in 

i i_ 

Eq. (|2"lj|) . However, that of (^,0^), which will be built from X \, Xfo, Ui, and Di, 
will contain more than one free parameter, and this will be incompatible with Eq.fl2.sj) 
as they have the same bulk S-matrix given by Eq.( |2.5[ ), Also interesting is to try to 
clarify the relation of the two classes of solution, distinguished by e = ±1, with the 
fused boundary S'-matrix of the kinks. 

We represent the boundary scattering |K a b(0i)K& c (02)) — ► |K C( ^(— 6' 1 )K^ e (— # 2 )) 
in Fig. 3. 




Fig. 3 Boudnary scattering of the two-kink states 



Let us begin with the bound states (ip^^n)- Using the fusion equations Eq.(^J), 
we can construct their boundary S'-matrix by combining Uq, Di and the bulk S-matrix 
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of the kinks as follows: 



4,^) = 2 (^0(01)00(02)5^(0! + 2 ) - C/ o (0i)^i(02)«S^(0i + 9 2 ) 
+D 1 (e 1 )D 1 (0 2 )S 1 £{9 1 + 9 2 ) - I>i(W(02)SiJ(0i + 02)) 

<*(*) = 2^(0o(0i)0o(0 2 )sg(0 1 + 2 )+ J D 1 (0 1 )[/o(0 2 )^J(0i + 2 
-^0(01)^1(02)5^(01 + 2 ) - -Di(0i)-Di (02)5^(01 + 2 ) 

2 1 2 1 



(3-4) 



= ^(0o(0i)0o(02)5p(0 1 + 2 ) + C/o(0i) J Di(02)5? 1 (0i + 2 
-£>i(0i)£/o(02)5!j(0i + 2 ) - -Di(0i)Z?i(0 2 )5i^(0i + 9 2 ) 



R% i4 ,{0) = | (^0(01)00(02)5^ (0i + 2 ) + C/ o (0i)A(02)5iJ(0i + 2 ) 

+D 1 (e 1 )u (e 2 )s 1 ^(e 1 + 2 ) + D 1 (e 1 )D 1 (0 2 )s\\(e 1 + e 2 )\ 

where the superscript 2 on the boundary S'-matrix refers to the second set of bound 
states. 

As an example the explicit computations for -R 2 ^,(0) and R 2 ^^ look as follows: 
TTJftATT n (ft n \ r 

Rip,ij(9) 



,2 , m _ 0o(0i)0o(0 2 )_ 

A0. 



a/2 (l + ,4sinh|-) (l + Asinhf 



< A0 
cosh - ( 1 - A 2 sinh 2 + A 2 sinh 2 



2 



« sinh — \ —1 — A 2 sinh 2 — + A 2 sinh 2 — 

2 V 4 2 



i? 2 = iaiAsiDheu {e 1 )u (e 2 ) ( co:Al A9 n , Ginh A ^ 

^ v^(l + ^sinh^) (l + Asinhf" 



Comparing them with the elements in the first column of the boundary S'-matrix 
given in Eq.(^l^), and after dividing by an overall factor 

aiAi7o(0i)0o(02) ( , A0 n . A9 n 

cosh — ■ i sinh — - — 



V2 (l + A sinh ^) (l + A sinh 



we see that the above RL ^ has the form 



cosh -G+(0)-i sinh -G_(0) (3.5) 
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with 



G±{6) = -(a 1 Ay 1 (cosh 



A9 r , 



i sinh • 



A9r 



±1 - A 2 sinh 2 ^ + A 2 sinh 2 



(3.6) 

Comparing Eq.( |3.6[ ) with Eq.( 2.1C| ), we deduce that this fused S-matrix corresponds 
to the e = 1 case since the coefficients of sinh 2 | are the same. Moreover, considering 
^-independent terms in both equations, we find 



a± = — {ol\A) 1 ( cosh 



A0 n 



i sinh ■ 



A9 n 



' / ±l-A 2 smh 2 ^ 



(3.7) 



Notice that the coefficient of the sinh 2 |-term in Eq.( |3.6| ) is given by 



Substituting the expressions for a± into this, we can produce the corresponding coef- 
ficient in Eq.( 2.10 ). In addition the a± satisfy 



a , 



a_ 



2 1 



sin 



n-y 
16 



which is consistent with Eq.( 2.1lj ). The same computations for the other two compo- 
nents 2? 2 0, -R^ ^ also lead to the same conclusion as above. We therefore confirm that 
this fused boundary S'-matrix indeed reproduces Eq.(2.9) with e = 1 and the mapping 
of the boundary parameter is given by 



%A 



a + + a_ 



l-(sin^)- 



(3.8) 



When (?/>!, 0^) scatter with the boundary, the out-states are linear combinations of 
,1 XI 



(tpn, 4>n) and two other states. For example we have 

A 



+ BC sinh 6 cosh ■ 



cosh i sinh — 

2 2 



|fii(-0)) 



1 



sinh 9 I cosh 



A6 



+ i sinh • 



A9 r 



.c 



B(l + s)+ i— (1 — s) (cosh 9 - cosh A9 r , 



\n 2 (-9)) , (3.9) 
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where ^, R\ ± are amplitudes to be given later. The two new states are given by 

|fii(0)) = -^(iKoi^Kio^-lK^^Ki,^))) , 

l«2(f)> = ^(|K o i(0 1 )Ki 1 (^)) + |K 1 i(^ 1 )Ki o (^ 2 ))) 

which are orthogonal to states. 

These new states do not scatter with (V'n'^n) a ^ an i n the bulk and their bulk S 1 - 
matrix is different from Eq.( [2.5D . Moreover, the boundary potential in Eq.( |2.8[ ) contains 
only the (t/j, 4>) pair and we should not expect any new particle to be created by the 
action of the boundary. Therefore, we should eliminate these extra states appearing in 
the fusion procedure so that the out-states are the linear combinations of only 
This is possible if we choose 

B = and s = 1 or C = and s = -1 . (3.10) 

The same conclusion is also arrived when the scattering of \ij)\{Q)) with the boundary 
is considered. 

With the above restriction we can reproduce Eq.([2.9|) by fusing the boundary and 
bulk S'-matrices of the kinks as follows: 




(3.11) 



Repeating the analysis given for ^ before on the above amplitudes for the two cases, 
we can show that these scattering amplitudes coincide with that obtained in Eq,(fO|). 
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We can summarize these in the following mappings: 



either B = s = l C= ^± — - 



sin ^)i-(sin^H 



(3.12) 

(sin^)t-( si n^)-t 



or C = s = -l B = ± 16 > ,\ 16 ' — e = l. 



The two classes (e = ±1) of the solutions presented in Eq.( |2.9[ ) are indeed compatible 
with the boundary S-matrices obtained from the bound states. 

Another aspect of the result is that the two N = 1 supermultiplets distinguished 
by the superscripts 1,2 scatter differently with the boundary since the scattering of 
(V^, <t>\) is given by the ^-matrices with e = ±1, while that of (V^, (ftn) by the S'-matrix 
with e = 1. Again we can require that the boundary potential should not add any new 
particle states in the theory and the identification of = {ip^} and |</>^) = |^) made 
in the bulk should hold for the boundary as well. This dictates that the boundary 
scattering amplitudes for these two multiplets be identical. Thus we can find that 



C = 0, s = 1, and iAB = * /sin — 

V 16 

with e = 1, and there is only one free parameter associated with the boundary. 



4 Boundary Supersymmetry 

The SSG model has a supersymmetry and the bound states (ip n , (j) n ) (we will drop the 
superscripts from now on) transform into each other under the action of the SUSY 
charges Q, Q in the following way |J: 

Q\M0)) = ^ 2 |i(S)), Q\M0)) = ^i^z~ 0/2 \MO)) , 

Q\M0)) = ^i^z d,2 \M0)) , QIM)) = v^e-^I^W), 

(4-1) 

where m n = 2 sin jq is the mass of the n-th bound state and we have not included the 
bound states of the SG sector since the SUSY charges act trivially on them. 

One can show that the bulk S'-matrix of the bound states is invariant under the 
action of the SUSY charges, namely, 

Si2(9)Q 12 (6) = Q 2 i(-e)S 12 (6) , S 12 (8)Q 12 (9) = Q 21 (-6)S l2 (9) . (4.2) 

It is interesting to examine the action of these SUSY charges on the boundary S-matrix 
to see whether SUSY can be maintained with the presence of the boundary. In fact, it 
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is argued that one can retain only 'half of the super symmetry Q ± Q in the presence 
of a boundary We shall see that this is the case for the SSG model. 

One can write the SUSY charges Q and Q as 2 x 2 matrices when acting on one- 
particle state as 

(9 in \ 
e2~~ \ 

8 | in I i 

e2+x 

(° 

V e 2 4 



Q 
Q 



_0_ \ ITT 

e 2^ 4 



(4.3) 



where we have arranged the basis in the order of ip n ,(f> n . 

Consider the action of the linear combination of charges Q(0) + (3Q{6) on the bound- 
ary S'-matrix: 

R{9) \q{9) + 0Q(6)] - \q(-6) + 0Q(-O)\ R{9) . 



Using Eq.(2.9), we deduce that for the above to vanish, we have either 

(3=1 for a + + q_ = 



(4.4) 



or 

P = -1 for a+ - a_ = . (4.5) 



From Eq. fl2.il ) we see that for the relation a+ = ±a_ to be true, we must have 



| ay. | — > co. In this limit, the boundary S'-matrix becomes diagonal 

/ cosh | + i/3sinh| 
R(e) = K (e)[ " I (4.6) 

\ cosh I - i/3sinh| 

and there is no boundary free parameter left, a conclusion which is in agreement with 
H that uses a different approach. 

We would like to comment that this SUSY preserving boundary S'-matrix, when 
regarded as the boundary reflection matrix of a lattice model, is the one that gives rise 
to a spin chain that possesses the quantum group Ui„2gl(l\l) symmetry. 



5 Discussion 



In this paper we considered the SSG model with boundary from two points of view. 
One is to consider the BYBE of the SSG breathers which are related to the eight- 
vertex free fermion model and the other is to use fusion of two kinks which are related 
to RSOS(4) model. By matching the boundary S-matrices obtained from these two 
approaches and requiring that there is only one N = 1 supermultiplet for the bound 
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states, we reduced the four parameters s, A, B, C in the SUSY sector of the SSG soliton 
scattering with the boundary to one, and also constrained the solutions of the BYBE 
to the free fermion model to that with e = 1. 

If we change to i(3 in the SSG model, we get the supersymmetric sinh-Gordon 
model with the difference that (?/>i , <pi ) are the only particles which are not soliton- 
antisoliton bound states as there are no solitons (anti-solitons) in the theory. Therefore 
the above argument for the sine-Gordon theory does not apply here and both of the e = 
±1 solutions are allowed. Since these solutions are mutually exclusive and associated 
with boundary potentials, they must be two different classes of the boundary potential 
which preserve integrability. In addition, each of these potentials should have at least 
one boundary coupling parameter which is related to that in the boundary S-matrix. 



The SSG boundary potential in Eq.2.8, derived from the condition that there exist 



the conserved charges at the first order, contains five parameters. If we compare this 
with the SSG boundary soliton S-matrix, which has only three parameters (two from 
the sine-Gordon ]l4| and one from the SUSY sector), we can think of two possibilities; 
One is that two of the five parameters in Eq.([2.8]) should disappear when higher order 
conserved charges are constructed. The other possibility is that all the five parameters 
survive and the soliton boundary S-matrix introduce two more parameters in its overall 
CDD factor. 

Furthermore, boundary SUSY is realized only when the five parameters are fixed 

as; 

A = ±8,M = ±l,0 o = O,e = e = O. (5.1) 

So there is no free parameters and the ± are related to the "half SUSY Q ± Q re- 
spectively. For the boundary 5-matrix, we found that indeed only "half SUSY can 
be realized. In which case, the boundary 5-matrix becomes diagonal given in Eq.(|4,6j) 
without any parameters, but it also comes in two classes which possess the "half SUSY 
respectively. This means there is no off-diagonal (0 — >• tp, t[? — > (f>) scattering amplitude 
and can be understood from the fact that the boundary potential deos not includes any 
fermion number violating term as far as e = e = 0. This SUSY, however, acts trivially 
on the sine-Gordon soliton sector, which contains two free parameters, this leads to the 
conclusion that at the SUSY points these two parameters become unphysical. 
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Note Added 



After we completed and submitted our paper, we have noticed that \22\ obtained a 
similar result as ours on the boundary supersymmetry covered in sect. 4. We thank the 
authors for the information. 
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Appendix: 

The Eight-Vertex Free Fermion Model with Boundary 



We present solutions to the boundary Yang-Baxter equation for the general eight-vertex 
free fermion model here. 

Recall the boundary Yang-Baxter equation takes the form [|l], || 

K 1 (u 1 )R 1 2(u 1 +u 2 )K 2 {u2)R2i{u2-ui) = Ri2(u2-ui)K2(u2)R 2 i(ui+u 2 )K 1 (u 1 ) (5.2) 

where R(u) and K(u) are, respectively, the bulk and boundary i?-matrices. 



The eight- vertex free fermion model has been studied by a number of authors [11 
15, [R], 17, 18, 19] the bulk iZ-matrix takes the form 



R(u) 



\d 



d \ 



(5.3) 



where a±,b±,c,d denote the usual vertex weights that depend on the spectral parameter 
u and other parameters of the model. These weights satisfy the free fermion condition 

a+a_ + - c 2 - d 2 = (5.4) 

and i?-matrices with the same parameters T and h given by 

2cd , a 2 _ + bl - a\ - b 2 _ 



h 



(5.5) 



commute. 

The extreme anisotropic limit of the bulk i?-matrix commutes with the quantum 
spin chain with local Hamiltonian given by 

h 
2 

Thus h can be interpreted as the bulk magnetic field of the spin chain. Using 
duality transformation or other argument, one can show that critical points must occur 
for \h\ = 1. Here we shall adopt the parameterization given in [15]; For \h\ < 1 the 
vertex weights are given by 



n 



3,3+1 



a 



(5.6) 



a± 
b± 
c 
d 



cosh 7 cn | ^ sinh 8 sn | dn | 
cosh S sn | dn | ± sinh 7 cn | 



dn| 

k sn I cn | 



(5.7) 
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with 

r = r /\ jr , k = tanh( 7 + <5) (5.8) 

cosh (7 + 0) 

where k is the modulus of the elliptic functions. While for \h\ > 1 R 



a± = cosh 7 dn | ± cosh 5 sn | cn | 

= ± sinh 7 dn | — k sinh <5 sn ^ cn ^ 

c = cn I dn I 

d = A;' sn I 



(5.9) 



with 



k' 

k sinh (7 + 5) 



/i = coth (7 + ,5) (5.10) 



where k is the complementary modulus of the elliptic functions. 

These i?-matrices satisfy the unitarity and cross-unitarity properties 

R(u)R(-u) oc 1 , R(u) tl R(2K - n)* 2 oc 1 . (5.11) 

So the crossing parameter is 2K a half-period of the elliptic functions. In general, 
due to the asymmetry, the i?-matrices do not have crossing symmetry. For the regime 
\h\ < 1, when 7 = 5, we have 

R(2K - u) t2 = afR(u)af , (5.12) 

hence the basis of the C 2 space can be given an interpretation of up-down spins. As 
for the regime \h\ > 1, when 7 = 0, we have instead 

R{2K - uf 1 = R(u) . (5.13) 

A natural interpretation of this property can be given if we regard the basis of the C 2 
space as two distinct degrees of freedom such as a boson and a fermion. 
Note also that for 7 = 0, these i?-matrices are regular in that 

R 12 (0) oc V 12 (5.14) 

where V12 is the exchange operator of the spaces 1, 2. 

With the bulk i?-matrix given, it is straight forward to solve for the boundary 



if-matrix using Eq^^). We again find that in order for solutions to exist, the bulk 
parameter 7 has to vanish, essentially this is due to the asymmetry of the vertex weights 
b±. With this restriction, the most general solution is given as follows: 

2 The spectral parameter here differs from that in jL5| by a shift of K. 
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For \h\<l, 

G + {u) cnf dnf -G_(«) sn f 2e sn f cn § dnf 



where 



2 sn | cnf dnf G+(«) cn f dn f + G_(u) sn f 

(5.15) 



^, /c ((1 — efccosh<5) ct + + sinh<5a_) 2 u 
G + - a+ + ecosh8-k sn 2 

k ((1 — ek cosh (5) a_ — k' 2 sinh (5a+ , 

<_r — Q!_ i r— r i Sn 77 

e cosh o — k 2 



Here a± are free parameters associated with the boundary that are related by 

k' 2 a 2 + + a 2 _ =2(V 1 cosh(5-e) (5.16) 

and e = ±1. 
For \h\>l, 

I G + (u) cnf +G_(n) sn f dn § 2e sn § cnf dnf 

K{u) = 

\ 2snfcnfdnf G+(u) cn f - G_(u) sn f dn f 

(5.17) 

where 

2 u a±e/c sinh 5 + a.^-k k cosh 5 
G ± («) = « ± - sn fc . + £fcsinh ^ • 

Again a± are free parameters associated with the boundary that are related by 

efcsinh(T 



and e = ±1. 



a 2 , - c?_ = 2e I 1 + 



fc' 



Hence for each regime there are two classes of solution distinguished by e = ±1 with 
one independent parameter. 

Below we shall list a number of properties of these if -matrices. 

1. Regularity 

K(0) oc 1 . (5.18) 
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2. Unitarity 

K(u)K{-u) oc 1 . (5.19) 

3. Corresponding Quantum Spin Chains 

With the boundary if-matrix and the bulk i?-matrix, one can construct an inte- 
grable quantum spin chain with fixed boundary condition using the prescription 
given in Q. The quantity that generates the spin chain Hamiltonian and other 
commuting conserved charges is 

t(u) = tr (k(K - u; a)T(u)K(u; a)T -1 (-u)) 

where T(u) is the usual bulk monodromy matrix constructed out of the i?-matrix, 
and a, a denote, respectively, the parameters a±,a± that are associated with the 
two boundaries of the spin chain. A new feature of the model is that tr (K{K; a)) 
vanishes, result we have 

t(u) = (const. )ul + (const. )u 2 7i + • • • , 

which is common in models with a supersymmetry.0 In the \h\ < 1 regime, we 
find 



U oc 



■ 7V-1 



cosh 2 5 — k 2 



/ cosh 5 + k 
cosh 5 — k 



+ 



' cosh 5 — k 



V V 



cosh 5 + k j j+1 



I sinh 5 + ik' 



I sinh 5 — ik' 



sinh 5 — ik 



7 Oa 



+ ^(a I +ea I] + 



v sinh 5 + ik' 3 
1 — ek cosh 5 — ik' sinh 5 ( a_ 



1 — ek cosh 5 H sinh 5 \ a+ 



ik \ a 



'N 



2(k + e cosh 5) 



a + (1 — ek cosh 5) + a_ sinh 5 
While for \h\ > 1, we find 



(ea N + (7+) 



(5.20) 



U oc 



-iV-l 



fc 2 cosh 5 -l 

3=1 



' k sinh 5 — k' 



' k cosh 5 + 1 
k cosh 5 — 1 



k sinh 5 + k 



' k cosh 5 — 1 
k cosh 5 + 1 



k sinh 5 + k' 



k sinh 5 — k' 3 3 



v v 



a 



3+1 



ft _i 



+ 



+ eaj 



+ 



+ l) ^JV + — (zo N + 4) 



(5.21) 



The same property is also observed in u.% for a special case of this model. 
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4. Quantum Group symmetry 

These spin chains include, as a special case, that obtained by [17] which has the 



quantum group Up 7 q(gl(l\l)) ( or CH q {2) ) symmetry where p, q are related to k, 5 



(see ]17], 18, T^]). This can be seen by taking the limit \a+\, |a_| |«-|) — > oo, 

which, from the relation satisfied by them, implies that the ratio of them is 
finite. For the regime \h\ < 1, taking 2p = ^ = ±ik gives the quantum group 

symmetric spin chain. While for \h\ > 1, we have to take = — ^ = ±1. Note 
that in this limit the i^-matrices are diagonal. 

5. Symmetric Limit 

It is well known that in the symmetric limit 5 = 0, ( in addition to 7 = ) the R- 
matrix for the regime \h\ < 1 is a special case of Baxter's symmetric eight-vertex 
model. Therefore, in the symmetric limit, the JC-matrix in the \h\ < 1 regime 
should also be a special case of that obtained in [|(], ^TJ. Recall that there are 
three boundary free parameters £, A, there, but we have only one free parameter 
here. We find that indeed our solution corresponds to the case where one of the 
above three parameters vanishes with the remaining two related. It is intriguing 
that this way of approaching the symmetric free fermion gives rise to K matrices 
of less boundary free parameters. 
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